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INTRODUCTION 

Aerospace vehicles are subjected to extremes in temperature and 
moisture (hygrothermal ) conditions within their flight regimes. These 
extremes in hygrothermal conditions have resulted in a rather intensive 
research effort on the part of materials engineers to understand their 
effects on the internal stress distribution and resulting behavior of 
resin matrix composites. Recent finite element investigations [1-3] 
have shown the significiance of a uniform temperature change on the 
interlaminar stress distribution in composite laminates with free 
edges. Pipes et al [4] performed a "unified" treatment of the hygro- 
thermal problem with a gradient rm>isture condition imposed; however, 
since lamination theory was employed the interlaminar effects were not 
assessed. 

The moisture diffusion work of Shen and Springer [5], using one- 
dimensional analysis, has shown the existence of transverse moisture 
concentration gradients in composite laminates. The works of Whitney 
[6] employing a three-dimensional analysis technique with average laminate 
properties and that of Farley and Herakovich [7] using a more general 
two-dimensional finite difference scheme have shown the possibility of 
large transverse (z) and lateral (y) moisture gradients, especially near 
free edge boundary layer regions. 

In the present study, a detailed investigation of the effect of 
uniform and gradient hygrothermal loading on the interlaminar stress 
distributions in finite width syimnetric composite laminates is presented. 
These stress distributions are compared with those produced by mechani- 
cal loads and the residual stresses induced during the curing process. 

Many cases were studied during the course of the investigation; however. 
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due to limitations of space only selected results will be presented to 
demonstrate the various phenomena. It will be shown that moisture 
induced stresses can be larger than the stresses due to thermal and 
mechanical loading, and that moisture gradients can significantly alter 
the stress distributions. 

The method of analysis used in this investigation was a computer 
implementation (NOLACS) of the displacement formulation of the finite 
element method. The program NOLACS is an improved version of NONCOM [8] 
and NONCOMl [9] which utilizes the full three-dimensional non-linear 
mechanical material properties including hygrothermal effects. A 
modular form of NOLACS utilizes low speed storage to greatly reduce the 
problem size restrictions of its predecessors. 

PROBLEM HPSCRIPTION 

The problem under consideration is the stress analysis of a long, 
finite width, symmetric composite laminate subjected to uniform axial 
mechanical loading and gradient hygrothermal loading (Fig. 1). Indi- 
vidual layers in the laminate are considered to be homogeneous ortho- 
tropic materials, and all stresses and strains are assumed independent 
of the axial (x) coordinate. The constitutive equation for such mater- 
ial behavior may be written in condensed i'sotation as 

{o} = [C]({e} - {a}T - {B)M) (1) 


where : 

[C] is the transformed 6x6 stiffness matrix 
{a} is a 6 X 1 stress vector 
{e} is a 6 x 1 strain vector 
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{a} Is a 6 X 1 vector of thermal coefficients 
fe) Is a 6 X 1 vector of moisture coefficients 
T is temperature 
M is moisture concentration 

The finite element formulation of the problem will not be presented 
here; interested readers may find the detailed presentation in refer- 
ences 8 and 9. Although previous applications of this formulation have 
been concerned only with the case of uniform temperature or moisture 
distributions, the initial assumptions in the problem formulation do not 
preclude its application to gradient conditions. 

Utilizing the formulation described above, several laminates have 
been examined using the finite element grid shown in Fig. 1 where 
symmetry has been invoked to limit the analysis to one-quarter of the 
laminate cross section. This model consists of 216 nodes with 376 
constant strain triangular finite elements. The synmetry displacement 
boundary conditions along z=y=0 are also shown in the figure as well as 
the stress free conditions along z=H and y=b. Although the stress free 
boundary conditions along the free edges will not be explicitly satis- 
fied due to problem formulation (displacement method), low number of 
degrees of freedom and the order of the finite elements used (constant 
stress elements), a sufficiently adequate stress distribution for 
relative comparisons and general conclucions is obtained. The laminates 
studied were [ 02 / 90215 , [90/02/90]^, [0/±45/90]^, and [90/0/±45]^. 

They were chosen because they exhibit significant interlaminar stresses 
in regions where hygrothermal gradients may radically influence these 
distributions. Both idealized one-dimensional uniform moisture gra- 
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dients (Fig. 2), and two-dimensional nonuniform gradients representing 
actual moisture distributions (Fig. 3) were applied to the laminates 
studied. 

The gradients representing actual moisture distributions were 
generated by a second order accurate two-dimensional hygrothermal finite 
difference program called HYDIP [7]. This finite difference program 
considers the problem of moisutre diffusion where the diffusion coef- 
ficient D may be a function of temperature T(y,z), moisture concen- 
tration M(y,z) and fiber orientation e, i.e. 

D = D(T(y.z). M(y.z). e) (2) 


The governing field equation 


+ D + + D ^ 

3y 3y y ^ 3Z 3Z z ^ 


iH 

3t 


(3) 


is solved by the finite difference procedure subject to the symmetry 
boundary conditions 


dz 


z*0 


= 0 ; ^ 
3y 


y*0 


and the variable surface conditions 


(4) 


M=M(z)|y,^; M=M(y)|^.„. (5, 

In generating moisture distributions the laminate was assumed initially 
dry with a uniform saturated (M = 1.4%) moisture condition imposed on 
the free surfaces. 

Although both programs NOLACS and HYDIP have the capability to 
allow tor hygrothermal dependent material properties and the finite 
element program has the additional capability of nonlinear material 
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behavior, the present study was limited to linear material behavior with 
the mechanical properties independent of hygrothermal conditions. All 
results were obtained using typical properties for graphite-polyimide. 
They are: 
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where the diffusion coefficients are given for 440'’K and the subscripts 
1, 2 and 3 correspond to the lamina material principal directions. 

RESULTS AND DISCUSSION 

Normalized stress profiles in the boundary layer region are pre- 
sented in Figs. 4-14. for a variety of laminates and loading conditions. 
The loading conditions have been coded for clarity as described in Table 
1 and all stresses have been "normalized" with respect to 1000 psi 
(6.895 MPa). 

uniform Hygrothermal and Meohaniaal Loading 

The results described in this section show comparisons between the 
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Interlaminar stresses In [ 02 / 902 ]^ and [90/0/±45]j laminates for uniform 
hygrothermal loading {M*1.4/6, T*280®F), extensional strain loading 
(ej^«1.0%) and a curing temperature (aT»-200®F). From Figs. 4-6 it is 
evident that elevated hygrothermal conditions can produce stresses which 
are larger in magnitude than those produced by other types of loading. 
For the laminates studied, the maximum interlaminar hygrothermal shear 
stress Ty^ is larger than the maximum interlaminar normal stress with 
the largest hygrothermal stresses being those of the bi-directional 
[02/902 ]j laminate (Fig. 4). As indicated in the figures, the results 
are a function of the interface in question as well as the material 
properties and the laminate stacking sequence. Scaling the results in 
Figs. 4-6 indicates that a uniform moisture concentration of approxi- 
mately 0.5% tends to cwnpensate for residual curing stresses associated 
with a 200®F temperature drop. Similarly, for the material properties 
used and the laminates examined, a moisture concentration of 1.0% is 
approximately equivalent to a temperature increase of 400® F. As these 
results show, moisture effects are a significant environmental con- 
sideration and the remainder of this paper will be directed toward the 
moisture problem. Thermal trends can be estimated from the moisture 
results since the problem formilations are the same. 

One-Dimeneional Uniform Gradient Moisture Loading 

A uniform temperature condition in a laminate can be reached in a 
reasonably short period of time because of the relatively high coef- 
ficient of thermal diffusivity; uniformly saturated moisture conditions 
can, however, take years to be attained at room temperature [5]. This 
leads one to believe that gradient moisture conditions are generally 
present even in controlled laboratory conditions. 


\ 

The results presented in Figs. 7-10 show interlaminar stress 
distributions for four moisture gradients which are uniform over the 
laminate. They are designated tY and ±Z and depicted in Fig. 2. The 
interlaminar stress distributions associated with uniform gradient 
moisture loadings are compared with uniform saturated moisture loading 
in the figures. The +Y and +Z moisutre gradients can be considered to 
simulate moisture absorption while the -Y and -Z gradients can represent 
desorption. By considering these Y and Z gradient conditions indi- 
vidually, one can more readily iMderstand the basic phenomena. 

The results shown in the figures indicate that the effects of the 
uniform gradient are localized in the boundary layer region with only a 
small shift in stress distribution from the uniform saturated condition. 
It is also apparent that absorption is more critical with respect to 
interlaminar stresses. This is because of the higher moisture concen- 
tration in the boundary layer region. Though not shown, it was noticed 
that during -Z gradient loading (desorption) a change in sign of the 
interlaminar shear stress occasionally occurred. These uniform gradient 
results are idealized cases which can, at best, yield trends that 
further the understanding of hygroscopic interlaminar stresses. As 
shown in references 6 and 7, actual moisture distributions have steep 
gradients in both the y and z directions. The effect of two-dimensional 
gradients are examined in the remainder of this paper. 


".\Jo-dimenP tonal Konunij'ovm C,yvidifi.nt Mnist'Ave 


uOC 


A series of moisture concentrations generated by HYDIP were applied 
to several laminates and interlaminar stress distributions were determined 
for each of three mositure concentration cases. The mositure contours 
are shown in Fig. 3 and the resulting stress distributions in Figs. 11- 
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14. Comparing these stress distributions to those previously examined 
and noting the location of the steep moisture gradients, It can be 
concluded that the stress distributions are significantly affected by 
the relative location of the steep gradients In the boundary layer 
region. After the steep gradients have passed through the boundary 
layer region, as In gradient case III, the Interlaminar stress distri- 
butions resemble those of the uniform moisture case. These results also 
indicate that the 'nterlamlnar stress distributions do not appreciably 
change for average moisture contents equal to or greater than eighty 
percent of the uniform saturated condition. This later observation may 
be especially important for experimental work since it Indicates that 
the Interlaminar stress distribution Is essentially fully established at 
an average moisture content which Is eighty percent of the saturated 
condition. 

It Is apparent from the figures that the Interlaminar normal stress 
distribution Is Influenced more by the location of the two-dimensional 
gradient than Is the shear stress distribution. The sinusoidal varia- 
tion of the normal stress for case II is unlike the distribution for any 
other loading condition previously Investigated. Such a distribution 
could have important implications for crack arrest should there be a 
positive negative variation in the sign of the noma! stress; such 
variations have been observed during this study. It is also apparent 
from the figures that the magnitude of the stress away from the free 
edge Is generally lower for case II, and that the slnuscxilal variation 
of the normal stress occurs when both y and z moisture gradients are 
present In the boundary layer. 
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CONCLUSIONS 


It has been shown that elevated hygrothermal conditions can be a 
significant factor in the determination of interlaminar stresses in the 
boundary layer region of composite laminates. Pertinent trends demon- 
strating the individual contributions of lateral and transverse gradients 
have been shown. By utilizing the diffusion program HYDIP, mositure 
distributions representing actual conditions have been analyzed. These 
results show that rather dramatic changes in the stress distributions 
may occur when two-dimensional gradients are present in the boundary 
layer region. After steep gradients have passed through the boundary 
layer (even though the entire laminate has not reached moisture equili- 
brium) the interalminar stresses are very similar to the uniform hygro- 
thermal loading condition. 

The results presented here show that hygrothermal as well as 
mechanical loading should be included in any stress analysis of com- 
posite laminates with free edges. Further, these results which have 
been determined for constant materials properties indicate that addi- 
tional study is needed in the area of hygrothermal de ident material 
properties and nonlinear material behavior in order to more accurately 
assess the load carrying capacity of composites. 
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Table 1 Stress Codes 


Superscript - denotes type of loading 
e - mechanical strain = 1.0% 

M - uniform moisture distribution M = 1.4% 

C - curing process aT = -200®F 
T - uniform temperature distribution aT = 280® F 
±Y, tZ - one-dimensional uniform moisture gradients 
I, II, III - two-dimensional moisture distributions 
Subscript - denotes interface as shown in sketch 
U - upper - first interface from top surface 
M - middle - second interface from top surface 
L - lower - third interface from top surface 
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FIG. i FINITE ELEMENT GRID WITH STRESS FREE AND DISPLACEMENT 
BOJNDARY CONDITIONS 
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FIG. 2 ONE-DIMENSIONAL UNIFORM MOISTURE 
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FIG. 7 INTERLAMINAR SHEAR STRESS PROFLES FOR 

UNIFORM MOISTURE GRADENTS IN A 0:>2/9O2] 5 
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FIG, 9 INTERLAMINAR SHEAR STRESS PROFILES FOR 

UNIFORM MOISTURE GRADIENTS IN A [90/02/90]s LAMINATE 




FIG. 10 WTERLA^MNAR NORM/SL STRESS PROFLES FOR 
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FIG II NTERLAMNAR STRESS PROFLES FOR TWO-DIMENSIONAL 
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INTERLAMINAR SHEAR STRESS PROFILES FOR 
TWO-DIMENSIONAL MOISTURE GRADIENTS IN A 
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